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Abstract. In a number of papers it was shown that there are one-dimensional 
systems such that they contain solutions with, so called, overcompressive sin- 
gular shock waves besides the usual elementary waves (shock and rarefaction 
ones as well as contact discontinuities). 

One can see their definition for a general 2x2 system with fluxes linear 
in one of dependent variables in This paper is devoted to examining 

their interactions with themselves and elementary waves. After a discussion 
of systems given in a general form, a complete analysis will be given for the 
ion-acoustic system given in |^. 

Keywords: conservation law systems, singular shock wave, interaction of sin- 
gularities, generalized functions 



1. Introduction 

Consider the system 

{gi{u)v + g2{u))t + (53(m)w + 9i{y'))x = 0. 



(1) 



where fi,gj, i — 2, ...,4, j — 1, ...,4 are polynomials with the maximal degree m, 
(m, w) = (u(x,t),v{x,t)) are unknown functions with a physical range Q, {x,t) e 
M X IR+. We shall fix the following notation for the rest of the paper: 



m m 



/.(y) = ^a,,fc/, .gj(2;) = ^6j,fe/, i = 2,3,4, j = 1,2,3,4. 

fc=0 fc=0 

There are cases when there is no classical solution to Riemann problem for the 
above system. Sometimes, there is a solution in the form of delta or singular shock 
wave. In |S] one can see when a system in evolution form (i.e. when /2 = w, 5i = 1 
and 52 = 0) permits a solution in the shape of singular shock wave. With the same 
type of reasoning and a more effort, one can give the answer to the same question 
in the case system (QJ. 

The aim of this paper is to investigate what happens during and after an in- 
teraction of a singular shock wave with another wave. After a general statement 
about new initial data taken at interaction point (of course, true for delta shock 
waves, too) in Section 3, we shall present a detailed investigation in the case of the 
system (so called ion-acoustic system) 

ut + {u^ - v)x = 

(2) 

vt + (mV3 - u), = 

given in 
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Definitions and concepts used here are from [2], based on the use of Colombeau 
generahzed functions defined in 11 . They will be briefly described in the Section 
2. If one is not familiar with these concepts, he/she can assume that a solution 
to the above system is given by nets of smooth functions with equality substituted 
by a distributional limit. The reason why the generalized functions are used is to 
give opportunity for extending the procedure in this paper for arbitrary initial data 
when a system posses singular or delta shock wave as a solution. 

Few interesting facts observed during the investigations of system ^ are arising 
a question about possibilities in a general case. Observed facts are: 

(1) The singular shock wave solution to a Riemann problem for (0) always has 
an increasing strength of the rate 0{t), t ^ oo. (The strength of the shock 
is a function which multiplies the delta function contained in a solution, 
s{t) in JZI))- After the interaction, the resulting singular shock wave is 
supported by a curve, not necessary straight line as before, and its strength 
can be an increasing, but also a constant or a decreasing function with the 
respect to the time variable. 

(2) When the resulting singular shock wave has a decreasing strength (this can 
occurs during an interaction of a admissible singular shock wave with a 
rarefaction wave), after some time it can decompose into two shock waves. 
This is a quite new phenomenon. 

The structure of this paper can be described in the following way. 

In the second section we will introduce necessary notation and give basic notions 
based on the papers ^^'^ IHl- 

In the third section, one can find a way how to continue a solution to the general 
case of system after an interaction point (Theorem^. The basic assumption is 
that a left-hand side of the first, and the right-hand side of the second wave can be 
connected by a new singular shock wave. The conditions for such a possibility are 
formulated trough a notion of second delta singular locus, see Definitional Explicit 
calculations for a geometric description of the locus are possible to perform for 
system but we shall omit it, to preserve readers attention on the further topics. 

The results given in these sections are used in the next one devoted to special 
case 

The first part of the fourth section is devoted to description of a situation which 
can occur after a singular shock and a shock wave interact. In the same way one 
can do the same for two singular shock waves, as one can at the end of this section. 

The final, 5th section, contains the most interesting and important results about 
singular shock and rarefaction wave interaction. In that case the decoupling of a 
singular shock into a pair of shock waves, already mentioned before, can occur. The 
analysis is done when a singular shock wave is on the left-hand side of a rarefaction 
wave. But one can easily see that these results can be obtained using the same 
procedure when a singular shock is on the other side of a rarefaction wave. 

2. Notation 

We shall briefly repeat some definitions of Colombeau algebra given in ^Jj and 
IH|. Denote M^^ := M x (0, oo), E^^ := R x [0, oo) and let C^{n) be the algebra of 
smooth functions on 17 bounded together with all their derivatives. Let C|?°(M^) be 
a set of ah functions u G C°°(M^) satisfying 'u|rx(o,t) G C^(R x (0,r)) for every 
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T > 0. Let us remark that every element of C^(K^) has a smooth extension up 
to the line {t = 0}, i.e. C^{M.%) = C^{^)- This is also true for C^(M^). 

Definition 1. fM,g(R+) is the set of all maps G : (0,1) x R, {e,x,t) l—^ 

Ge{x,t), where for every e e (0, 1), Ge G satisfies: 
For every (a, /3) e Nq and T > 0, there exists N eN such that 

sup \d^d^G^{x,t)\ = 0{e-^), as £ ^ 0. 

(x,t)eMx(0,T) 

^M,g(IR+) is an multiplicative differential algebra, i.e. a ring of functions with the 
usual operations of addition and multiplication, and differentiation which satisfies 
Leibniz rule. 

J^giM-X) is the set of all G € SM,g0^l), satisfying: 
For every (a, /3) e N^, a e M and T > 

sup \d^dfGe{x,t)\ - Oie"), as e 0. 

(x,t)eMx(0,T) 

□ 

Clearly, A/'j,(R^) is an ideal of the multiplicative differential algebra £M,g{^\), 
i.e. if Ge e Afg{Rl) and H, e 5M,g(K+), then G^He G A/'g(ffi+). 

Definition 2. The multiplicative differential algebra Gg{R\) of generalized func- 
tions is defined by GgiM.^) = £M,giR\)/^fg{R\)■ All operations in GgiM.^) are 
defined by the corresponding ones in £M,g{^+)- O 

If C^(R) is used instead of C^(M^) (i.e. drop the dependence on the t vari- 
able), then one obtains £M,g(R), ■^gi^), and consequently, the space of generalized 
functions on a real line, (^^(M). 

In the sequel, G denotes an element (equivalence class) in Qg{0,) defined by its 
representative G^ G £M,g{^)- 

Since C2°(IR^) = C^(R^), one can define a restriction of a generahzed function 
to {t = 0} in the following way. 

For given G G CJg(IR^), its restriction G\t=o G ^g(IS) is the class determined by 
a function Ge{x,0) G £M,g{^)- In the same way as above, G{x — ct) G QgiR) is 
defined by Ge{x - ct) G £M,g{^)- 

If G E Qg and / G C°° (R) is polynomially bounded together with all its deriva- 
tives, then one can easily show that the composition /(G), defined by a represen- 
tative /(Ge), G G Gg makes sense. It means that f{Ge) G £M,g if G^ G £M,g, and 
f{G,)- f{H,)eNgiiG,-H,eNg. 

The equality in the space of the generalized functions Gg is to strong for our 
purpose, so we need to define a weaker relation called association. 

Definition 3. A generalized function G G Gg{^) is said to be associated with 

u G 2?' (17), G ~ M, if for some (and hence every) representative G^ of G, G^ ^ u 
in 'D'{Q) as £ ^ 0. Two generalized functions G and H are said to be associated, 
G « iJ, if G — i? « 0. The rate of convergence in V with respect to s is called the 
rate of association. □ 

A generalized function G is said to be of a bounded type if 
sup |Ge(x,t)| = e»(l) as £ ^ 0, 

(x,t)eMx(0,T) 
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for every T > 0. 

Let u e 2?^oc(M). Let An be the set of all functions (j) e Cq°{R) satisfying 
^{x) > 0, X e R, / (j){x)dx = 1 and supp^ C [—1, 1], i.e. 

Ao = {(f>eC^ : (Va; e R)(t){x) >0, j (f>{x)dx = 1, supp^!> C [-1, 1]}. 

Let = e-V(^/e), xeR. Then 

where u*(f)e/Mg denotes the equivalence class with respect to the ideal Ng, defines 
a mapping of ■D^oo(lR) into Qgi^), where * denotes the usual convolution in V . It 
is clear that commutes with the derivation, i.e. 

Definition 4. (a) G € Qg(R) is said to be a generalized step function with 
value (t/o, yi) if it is of bounded type and 

GM - 

Denote [G] := yi — yo- 

(b) D G Gg{R) is said to be generalized split delta function {S6 -function, for 
short) with value {ao, ai) if D = aoD~ + aiD~^, where ao + ai = 1 and 

DG « {yoao + yiai)6, (3) 

for every generalized step function G with value {yo, yi). 

(c) Let m be an odd positive integer. A generalized fmiction d € Qg{R) is 
said to be m' -singular delta function {m' SD-function, for short) with value 
(/3o,/3i) iid = M- + (3id+, (3^-^ + = 1, rf± e ^^(M), {d^f « 0, 
i e {1, . . . , m - 2, m}, (rf^)™"! w 5, and 

rf™-iG«(rf-i+yi/3r')'^, (4) 

for every generalized step function G with value (j/o, 2/1 )• 

(d) Let m be an odd positive integer. A generalized function d £ Gg{R) is 
said to be m-singular delta function {mSD- function, for short) with value 
(/3o,/3i) if = M- + I3id+, 13^ +I3r = 1, d^ e GgiR), {d^Y « 0, 
i e {1, . . . , m - 1}, (d±)™ w (5, and 

rf'"G«(t/o/3o"+yi/3r)^, (5) 
for every generalized step function G with value {yo,yi)- 

□ 

In this paper we shall assume the compatibility condition Dd ~ 0, where D is 
S5- and d is mSD- or m'Sd-function. 

Suppose that the initial data are given by 

I \ uq, x < X 1 ^u, u- \ ^j. 

u\t=T = { ^ v\t=T = { ^ (6) 
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Definition 5. Singular shock wave (DSSW for short) is an associated solution to 
(PJ with the initial data © of the form 

u{{x - X), {t - T)) = G{{x -X)- c{t - T)) 

+ S{t){aQd-{{x - X) - c{t - T)) + aid+{{x ~ X) - c{t - T))) 

v{ix - X), {t - T)) = H{{x - X) - c{t ~ T)) (7) 

+ s{t){(3oD-{{x -X)- c{t - T)) + f3iD+{{x - X) - c{t - T))) 

+ ~s{t){^^d- {{x -X)- c{t - T)) + jid+iix -X)- c{t - T))) 

where 

(i) c G M is the speed of the wave, 

(ii) s{t), s{t) and s are smooth functions for i > 0, and equal zero at t — T. 

(iii) G and H are generalized step functions with values {uq,ui) and {vq,vi) 
respectively, 

(iv) di = aod^ + Pid^ and d2 = ^od^ + 71 are mSD- or m'SD-functions, 

(v) D = af)D~ + aiD^ is an S(5-function compatible with d. 

The singular part of the wave is 

s{t){at)d'' + ai(i+) 
_s(t){l3oD~ + l3iD+) + l(0(7od" + 7id+)_ ' 

The wave is overcompressive if its speed is less or equal to the left- and greater or 
equal to the right-hand side characteristics i.e. 

A2(wo,i'o) > Ai(mo,i'o) > c > A2(fii,'i;i) > Ai(ui,'i;i). 

□ 

Remark 1. (a) In [Sj one can find special choice for S5- and and d is mSD- or 
m'Sd-functions. For example S are given by the representatives 

mSD- and m'SD-functions can be chosen in the same manner. 

(b) Compatibility condition for an S(5-function D and an mSD- or m'SD-function 
d is automatically fulfilled if 

supp df n supp — supp d^ n supp = 

(c) Idea behind the above definition of products (0), Q and Q is the following. 
Starting point is that we know nothing about infinitesimal values of the initial data 
(carried on by step functions G and H above) around zero, but only that any such 
unmeasurable influence stops at the points ±e. The above mentioned definitions 
are made in order to get uniqueness of all products where step functions, Sd-, 
mSD- and m'SD-functions appear. With an additional information for and Hg 
around zero, one can choose D and d much more freely. For example, in G^ and 

are monotone functions (which is quite natural assumption), relation Q can 
be substituted by 

DG « 7(5, 7 can be any real between min{yo, yi} and max{yo, yi}- 

The possibilities in Colombeau algebra are even wider for specific systems instead 
of general case (0). One can look in for a good review of such possibilities. We 
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dealing with a system in a general form and it is the reason for using the above 
definition. 

(d) Due to absence of known additional facts for the general case ^ (hyperbol- 
icity, additional conservation laws,...), one can use the over compressibility as an 
admissibility condition. 

Definition 6. The set of all points {ui,vi) G such that there exists an singular 
shock wave solution (called corresponding DSSW) to Cauchy problem (|llti|) is called 
delta singular locus. We shall write {ui,vi) £ T)SIj{uq , vq) . If the corresponding 
DSSW is overcompressive, then it is called overcompressive delta singular locus. We 
shall write {ui,vi) £ DSL*(uo,fo)- D 

In the sequel, the term "solution" will denote generalized function which solves 
a system in the association sense. 

3. The new initial data 

Suppose that system posses a DSSW solution for some initial data. Assume 
one of the following. 

(i) If an mSD-function is contained in the above DSSW, then assume 

deg(5i) < m - 1, dcg(.g2) < m, deg(/2) < m. (8) 

(ii) If an m'SD-function is contained in the above DSSW, then assume 

deg(ffi) < m - 2, deg(g2) < m - 1, deg(/2) < to - 1. (9) 

Take the new initial data 

juQ, x<X Jwo, x<X 

\ui, X > X \vi, X > X 

for system (^, where ^ is a non-zero real. 

Definition 7. The set of all points {ui,vi) £ fl such that there exists an DSSW 
solution (called corresponding DSSW) to Cauchy problem (|1I10|I for some ( is called 
second delta singular locus of initial strength for [uq, vq). We shall write (ui, vi) £ 
SDSL^(uo, i^o) If the the corresponding DSSW is overcompressive, then it is called 
overcompressive second delta singular locus, and write {ui,vi) £ ST)Sh'^{uo,vo). □ 

Before the main theorem, let us give a useful lemma. 

Lemma 1. Suppose that (ui,vi) £ DSL(uoii^o)- Then (ui,vi) £ SDSL^(uo, wo), */ 

C>o. 

// the corresponding DSSW contains mSD-function, and m is an odd number, 
then the statement holds true for every real C. 

Additionally, (3i, i = 1,2, from Definition\^for the corresponding DSSW do not 
depend on ^. 

Proof. We shall give the proof for a DSSW containing mSD-function . The other 
case can be proved in the same way. 
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Inserting functions u and v from Q into system with initial data and 
taking account relations (|H1) or jnj, one gets 

f2{u) «/2(G) 

9i{u) ~9i(G) 
92{u) «52(G) 

/3(") ~h{G) + S(t)"-i(uia™-id- + uoa5""'d+)ma3,™-i 
+ + <d+)a3,™S(t)"a3,™^ 

g3(u) «53(G) + s(t)""'(^^i/3o""'d" + uo/3r-M+)m63,™-i 

g4(u) «/4(G) + s(t)'"(/3™-id+ + /3r"'rf+)64,™S(<)'"64,™<5 

There are two possible cases. Either s ^ and a^.m = &3.m = (i.e. deg(/3) < m—1 
and deg((73) < to — 1), or s = 0. In both the cases, the procedure which follows is 
the same, so take s ^ for definiteness. From the first equation of ^ one gets 

{f2{u))t + {h{u)v + f^{u))^ 

^-c{[h{G)] + [h{G)H + U{G)])5 

+ S(<)'"-i|(t)(7.i<-So + uo<"'7i)ma3,™-i<5' 
+ (/3(wo)/3o + h{ui)Pi)5' + S{tr6' « 0. 
One immediately gets the speed of DSSW, 

_ [/3(G)g + /4(G)] 
[/2(G)] 

and the relations 

Kis{t) = S{t)"'-^S{t) and K2s{t) = s(i)™, 
for some reals ki and K2. Finally, one gets 

ki(uiq;™"Vo + ■tioQ;™"Vi)"^a3,m-i + /3(uo)/5o + h{ui)fii + ^264,™ = 0. (11) 
Inserting all these relations into the second equation, one gets 
{gi{u)v + g2(u))t + {g3{u)v + g4u))x 
^{^cMG)H + .g2(G)] + MG)H + g,{G)] + s'{t)ig,{uoPo + gi{ui)Pi))S 
+ s{t){gi{uo)(3o + gi{ui)(3i + 53(wo)/3o + 53('"i)/3i 

+ Ki(uiQ;™"Vo + UoQ;™"Sl)'7^63,m-l + K2fe4,m)<5' ~ 0. 

The function s must be a linear one, say s'(t) — a, and the above functional 
equation gives the last two equations in M, 

-c[gi{G)H + 52(G)] + b(G)i/ + g4(G)] + a{gi{uo)Po + 5i(ui)/3i) = (12) 

and 

- c((5i(uo) + g3{uo)f3o + (5i(ui) +53(mi)/3i) 

(13) 

+ ki(uiq;™ 7o + uoar 7i)'7i63,m-i + ^264,™ = 0. 

In the above equations, only important fact about s is its derivative. Thus one 
can safely put s{t) = at + ( and if the above system pill3|l has a solution, then 
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{ui,vi) € SDSLf (mo, t'o) provided that s and s can be recovered. This is certainly 
the case when C > 0. If m is an odd number, then s — s(t)^/™ and s = s are always 
determined. 

The second part of the assertion, that f3i, i = 1,2 are independent of C is obvious 
from the above. □ 

Remark 2. From the proof of the lemma one can see that it is actually possible for 
C to take negative values, i.e. it is enough that C > — s(T), where T is a time of 
interaction when new initial data are given. 

The following assertion is crucial for the construction of weak solution (a solution 
in an associated sense) to after an interaction: At an interaction point of a 
DSSW and some other wave one can consider the new initial value problem which 
contains delta function. 

Suppose that the initial data are given by 

{uq, X < a I wqj X < a 

Ml, a <x <b and 0) — Ivi, a < x <b (14) 
U2, X > b [w2, X > b 

such that there exist a singular shock wave starting from the point x — a and a 
shock wave (or another singular shock wave) starting from the point x ~ b, a < b. 
They can interact if ci > C2, where Ci is the speed of the i-th wave, i — 1,2. For 
the simplicity we shall assume that 6 = 0. 

Let (X, T) be the interaction point of the overcompressive singular shock wave 
starting at the point x = a 

v}{x,t) =G^(x — cit — a) + (t) (^a^d^ {x — cit — a) + a\d'^{x — cit — a)j 
v^{x,t) =H^{x - cit -a) + s^{t)(^f3lD-{x - at - a) + (3lD+{x - cit - a)) (15) 

+ (t) (lod'' ix-cit-a)+ 7i d+ {x - Cit - a)j 
and the admissible (singular) shock wave 

u^(x,t) ^G^{x ~ C2t) + d^{t){ald^[x - C2<) + a?d+(a; - csi)) 

V^{x,t) ^H^{X ~ C2t) + s'{t)[plD-{x ~ C2t) + PlD+{x - C2t)) (16) 

+ s^{t){^ld^ {x - C2t - a) + 7i(i+(a; - C2t - a)j 

where G^, G^, and are the generalized step functions with values (mo,ui), 
(ui,M2), {vo,vi) and (wi,W2), respectively. Also, {al)"^^ + {a\)"^^ = (70)"' + 
(^*)™i = /3q + /3]; = 1, i = 1, 2. Here, mi = m if singular part of singular shock 
wave is mSD-function and mi = m — 1 in the case of m'SD-function. If the second 
wave is a shock one, then one can put = = = 0. 

The speed of a singular shock wave (as well as for a shock wave) can be found 
using the first equation in because of assumptions or jnj. For the first 
singular shock wave (|15|l we have 

{f2{u))t + (/3(")^^ + fA{u))^ « {h{G))t + {h{G)H + /4(G)), + (const si(i)<5), 
«(-ci[/2(G)] + [h{G)H + J^{G)])6 + const si(t)<5' « 0, 
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where the term const s^(t) is determined, but wc shall not write the exact value 
since it is not needed for the assertion. Missing argument in the above expression 
is a; — cit — a. 

Let Ti — {x = cit + a} and r2 = {a; = C2t\. Then [-Jr; denotes the jump at the 
curve Ti, z = 1, 2. Thus, one can see that the speed of that singular shock wave has 
the same value as in the case of shock wave, 

^ [h{G)H + /4(G)]r. 
[/2(G)]r, 

Also, 

[h{G)H + /4(G)]r. 



C2 



[/2(G)]r. 

Finally, one can see that the waves given by H15|l and H16() will interact at the 
point (X, T) if a < and ci > C2 , where 

a[/2(G)]r,[/2(G)]r, 



T = 
X = 



[/3(G)i/ + /4(G)]r,[/2(G)]r, - [/3(G)i? + /4(G)]rJ/2(G)]r, 

a[/3(G)g + /4(G)]r,[/2(G)]r, 

[/3(G)7J + /4(G)]r,[/2(G)]r, - [/3(G)i/ + /4(G)]rJ/2(G)]r. 



Denote by {u{x, t), v{x, t)) a solution before interaction time t — T consisting of 
waves H15I16|I . 

Remark 3. In the case of system one can easily calculate speeds of the above 
shocks and coordinates of the interaction point. The speeds of singular shock and 
entropy shock wave are 

u\- Vl-ul + Vo U^- V2- u\+ Vl 

Cl — and C2 = . 

Ml - Mo 1*2 - Ui 

If Cl > C2, then one gets 

X = ^^andT=^— . 



C2 - Cl C2 - Cl 

for the interaction point (X, T). 

Theorem 1. Let system 0) be given. Suppose that {u2,V2) G SDSL^(uoi ■f^o); 
C ~ (^1 + C2)/(3i(j^o)/3o + 3i(ui)/3i), where the constants i — 1,2, are defined by 

gi{u^)v'^ + g2{u'^)\(t=T) ~ Ci^(x,T) 
gi{u^)v^ + g2{u^)\(t=T) ~ C2<5(x,T)- 
The corresponding DSSW, {u,v){x,t) is given by 
u{x,t) =G{x - X - c(t ~ T)) 

+ S{t)(aod-{x -X -c{t~ T)) + aid+{x - X - c{t ~ T))) 
v{x,t)=H{x-X-c{t-T)) (-^7^ 
+ s{t){j3QD-{x -X-c{t- T)) + l3iD+{x - X - c{t ~ T))) 

+ l{t){jod-{x -X-c{t- T)) + -fid+{x-X - c(t - T))) 
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for t > T. By Lemma^ /3o and Pi are determined independently on so the 
definition of DSSW makes sense. 

Then there exist a solution to in the association sense such that it equals 

{u, v){x, t) for t < T — e, and it equals (u, v){x, t) for t > T + e. 

Proof. Take a constant to such that singular parts of the waves {ul{x,t),vl{x,t)) 



and {uj{x,t),v^{x,t)) has disjoint supports (i.e. cit- 



C2t > 4£, for t < T-toe, 



if one uses the construction of the SS, mSD and m'SD- functions defined above). 
Let us denote 





= {(x,t) : 


\x- 


X 


<toe 


+ e 


\t 


-T\<toe + e}, 


Ae 


= {(x,t) : 


\x- 


X 


< he, 


\t- 


-T 


< he}. 


Ae 


^{{x,t) : 


\x- 


x\ 


< toe- 




t = 


--T-toe-e}, 


Be 


^{{x,t) : 


X — 


X 


f <oe + 


e. 


t 


T\ <toe + e}. 


Ce 


= {{x,t) : 


\x- 


x\ 


<toe- 




t = 


--T + tae + e}, 


De 


= {(x,t) : 


X — 


X 


-toe- 






T\ <toe + e}. 



Define a cut-off function ^e(x,t) which equals zero for {x,t) S A^ and 1 for 
{x,t) G Ae. Let 



{utemp, Vte'mp){x, t^ 



{u{x, t),v{x, t)), t <T 
{u{x, t),v{x, t)), t > T. 

We shall prove that the generalized functions u and v represented by 

Ue{x,t) = Utemp{x,t)^e{x,t), aiid Ve{x,t) = Vtenip{x,t)S,e{x,t), X G R, t>Q (18) 



solve in the association sense. 
Denote 



Y{u,v) = 



h{u) 
gi{u)v + g2{u) 



and G(u, v) 



f3{u)v + fi{u) 
gz{u)v + gi{u) 



We have 



/ / F(u, v)^t + G(u, v)-ii^dxdt 
J Jr\ 

-. jj r(u, w)*t + G(?i, v)^xdxdt 



1p2 



F{u, u)*t + G{u, v)^xdxdt. 



for every test function * = G C(f (M^). 
The measure of the set is ©(e^), as e ^ 
|jF(u,«)*t + G(u,w)*,|ii..(R 
due to the assumptions in Definition 0] Thus, 

jj Y{u,v)-^t + G{u,v)-^^dxdt 



0, while 

M < const £-1+1/" 



-l/m 



0, as e ^ 0. 
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Using the divergence theorem for the second integral one gets 
/ / F(m, v)^t + G{u, v)^xdxdt 

= / F(u,w)*j^t+ j G{u,v)-^i'xds 
JdA^ J 

-II ¥{u,v)t'^ + G{u,v)x'i!dxdt. 

The last integral in the above expression tend to zero as £ ^ since {u, v) solves 
in \ Ag due to the construction. For the other integral one gets 

/ F(u,w)*i't+ / G(m, u)^'z/^ds 
JdA^ J 

= / F{u,v)'i'dx- / F{u,v)'i'dx+ / G{u,v)-^dt- / G{u,v)'ifdt. 

J Jb^ 

Functions and are i°°-bounded uniformly in e on the sides and D^. Since 
their lengths are ©(e), integrals over them tends to zero as e — > 0. 
Using the fact that f2{de) ~ one gets 



lim F{u, v)\t=T 

e— >0 







_(Cl + C2)S(x,T)_ 

as well as the construction of SS- and m'SD (or mSD)-functions, one gets 



lim / F{u^,Ve)dx 



Thus, there has to be true that 



lim / F{u^,Ve)dx 





C1+C2 





C1+C2 



*(X,T). 



This implies f2{u)\(x,T) ~ ^ 

9l{u)v + 92{u)\{X,T) ~ (Cl +C2)'5(X,T)- 



(19) 



Due to conditions © or @ one immediately gets /2(w)|(xt) ~ 0. Put ^ = (^^ + 
C2)/(5iK)/3o+3i("i)/3i). Then 

gi{u)v + g2{u)\t=T ^GH + G + s{T){gi{uo)/3o + gi{ui)(3i)5{X) 
and after another restriction on a; = X, 

gi{u)v + g2{u)\{X,T) ~ (Cl + C2)'^(X,T)- 

This concludes the proof. □ 
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Remark 4. The distributional limit of the result of the interaction is given by 

uq, X < cit — a, t < t 
ui, cit ~ a < X < C2i, t < T 

u{x,t) = < U2,X> C2t, t <T 

Uq, X < Ct + X, t > T 

U2, X > ct + X, t > T 
' Vq, X < Cit — a, t < t > 
Wi, Cit — a < X < C2t, t < T 

V{x,t) = < V2,X > C2t, t <T 

Vo, X < ct + X, t>T 
, V2, X > ct + X, t > T 

where Si = {(2^,^) : x = cit + a, t G [0,T], S2 — {ix,t) : x — C2t, t e [0,T] and 
S = {{x,t) : X — X = c{t — T), t E [T, 00). If the second wave (I16|) is a shock one, 
then S2 = 0. 

The above solution is continuous in t with values in This fact can be 

used in the approach similar to where the variable t is treated separately, i.e. 
when system ^ is considered to be in evolution form. 

The theorem shows that after an interaction of a singular shock with some shock 
or another singular shock the problem reduces to solving system with the new 
initial data 

Remark 5. (i) The solution to the interaction problem from Theorem^ is always 
associated with a lower association rate than the solution of the original Riemann 
problem. For specific system it seems possible to make more sophisticated con- 
struction in order to improve the rate. 

(ii) It appears that df are unavoidable correction factors even their distributional 
limit equals zero. 

The conditions ^ and ^ ensures that the new initial data at intersection point 
do not depend on mSD- or m'SD-functions in the solution. We have used them 
because the real nature of mSD- and m'SD-functions is not so clear yet. 

The above theorem will be used in the rest of the paper for investigation of 
interactions between singular shock waves and other types of waves in the special 
case of system ^ . 

4. Applications 

Consider now system ^ which a special case to The authors of IB] defined 
and proved existence of singular shock wave solutions for some Riemann problems 
of this system. 

In the present paper, we will investigate interactions of such solutions with the 
other solutions to Riemann problem for 10). In order to familiarize a reader with 
the presented results, let us give some basic remarks about such solutions. 

For a given Riemann data (uq^vq), {vo,vi), there are three basic solution types: 
(a) Shock waves 



i{x,y) = 



Uq, X < ct 
Ui, X > ct 



v{x,y) = 



Vo, X < ct 
Vi, X > ct 



(20) 
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where c ~ [u^ ~ ''^l/M and [ui^vi) lies in an admissible part of Hugoniot 
locus of the point (uq, wq)- 

(b) Centered rarefaction waves 

Uq, X < {uq — l)t 
u{x, t) = { x/t + 1, {uq ~l)t<x< (ui - l)t 
Ui, X > (ui — l)t 

(21) 

vo, X < (uo- l)t 
v{x, t) = { {x/tfl2 + 2xlt + Ci, (uo ~l)t<x< (ui - l)t 

Vl, X > {ui — l)t 

(1-rarefaction wave), where Ci = vq — Uq/2 ~ uq ~ 1/2, when (mi,i'i) lies 
in an 1-rarefaction curve starting at the point (uq^vq). Or 

Uq, X < (itQ + l)t 
u{x, t) ^ { x/t - 1, {Uo + l)t<X< {Ui + l)t 
Ui, X > (ui + l)t 

(22) 

Vq, X < {uq + l)t 

v{x, t) ^ { {x/tfl2 - 2xlt + C2, (uo + l)t<x< (ui + l)t 

Wl, X > {ui + \)t 

(2-rarefaction wave), where G2 — vq — uf^/2 + uq — 1/2, when {ui,vi) lies 
in an 2-rarefaction curve starting at the point (uo,wo)- 

(c) Singular shock waves (see Definitional) of 3'SD-type, 

r Mo, x < ct^ 

u{x, y) = x> ct\^ s{t){(^od, {x - ct) + aidj{x - ct)) 

: (23) 

(VQ,X<ct^ 

v{x,y) 2: > ct} + *(*)(/5oi57(a; - ct) + PiD^ix - ct)), 

where c — [u^ — v]/[u], and all other terms are as in that definition. That 
means 

D, « 6, (df)' « 0, i = 1, 3, (d±)2 « 6, (24) 

while {ui,vi) lies in a region denoted by Qt in 6 of the point (uo,fo) (see 
Figure 1). 

For an arbitrary Riemann problem to one can construct a solutions by the 
means of these waves or their combinations ([H])- 

While interactions of the first two types can be handled in a usual way, interac- 
tions involving singular shock waves are quite different and far more interesting, so 
they become a topic of this paper. 

The procedure for the singular shock wave interactions can be also used for sys- 
tems But a complete after-interaction solution highly depends on a particular 
system. That is the reason why we treat system (j^J only. 

In order to simplify notation, we shall substitute the point {X, T) in H10() by 
(0, 0) and then solve the Cauchy problem l|2ll()|l . 

There are no multiplication of v with u in system so in the sequel it will be 
enough to take ~ ao{t) :— aos{t), ai{t) := ais{t) and P{t) :— s{t), i.e. to 
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look for a solution of the form 

u — G{x — ct) + {aci{t)d^ [x — ct) + ai{t)d^{x — ct)) 
v^H{x-ct)+(3{t)D{x-ct), 

where G and H are generalized step functions, while d is 3'SD- and D is S(5-function 
and c e M. 

Let us determine SDSL of (0) for some {uq, vq) e M^. 
Substitution of H25|l into the first equation of the system gives 

_ ul - Vi - ul + Vq 

Ui ~ Uq (26) 

where c is the speed of the wave. After neglecting all terms converging to zero as 
e — > 0, the second equation becomes 

dtH,{x - ct) + (3'{t)5{x - ct) ~ cl3{t)6'{x - ct) + d.^{\Gl - G) 
+ {uial{t) + uoal{t))S'{x - ct) = 0. 
Thus, the following relations has to hold. 

(3'{t) = c{vi - vo) - {^ul -uo- ^ul + Ml) =: fc, (27) 

i.e. 

(3{t) = kt + (, since /3(0) = C 

and 

uial {t) + unaj (t) = cP{t) . (28) 
Like in 1^ one can see that the overcompressibility means 

Uo — 1 > c > ui + 1, 

i.e., Vi lies between the curves 

D = {{u,v) : V = Vq + + U — UqU — Uq} 
E = {{u, v) : V = Vq — U + UqU — Uq + Uq} , 

and Uq — ui > 2. 

Denote by Ji the union of the parts of admissible Hugoniot locus 

o (f \ ( .fUQ + u [ (uo - w)^\\ 
S'i = |(u,u): v-vq^[u~uq)\^ — - — + yi — j|, 

and 

Q u \ I -.f^o + u [ {uq -u)'^\\ 
{[ui,vi) : v-vq = [u-uq)\^ — yi — j|, 

for u €[uq — VT2, Uo — 3] . Note that Si is not an ith shock curve but only a label. 

The points between the curves D and E, and on the left-hand side of Ji defines 
the area denoted by Q7 in Here, this area is called delta singular locus. 

One can easily check that system (|26I28|) has a solution if and only if f3{t) > 0. 
Depending on fc, defined in H27|) . there are three possibilities for a resulting wave: 

(i) If fc > 0, then /3'(i) > and {ui,vi) G Qj. The resulting singular shock has the 
same properties as before, i.e. its strength increases with the time. 
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Figure 1. 

(ii) If fc = 0, then (3 = const — > Q and the corresponding part of a singular 
overcompressive locus is Ji . The result of the interaction is a new kind of singular 
shock wave, its strength is a constant with respect to the time. 

(iii) If A: < (this means that the point {ui,vi) is on the left-hand side of Ji), 
then the resulting singular shock wave has much more differences from the usual 
one (with an increasing strength). Its initial strength equals Cj /3(0) = C > 0, but 
linearly decreases in time. At some point Tq the strength of the singular shock 
equals zero and the singular shock wave does not exist after that. In the rest of the 
paper we shall see some cases when this happens. The new initial data for time 
t = Tq are the Ricmann ones, and the solution after that time can be find in the 
usual way, by using the results in (6j. 

All the above facts are collected in the following theorem. 

Theorem 2. The SDSLq, C > 0? for I^S110\} is the area bounded by the curves D, E, 
S2 \ Ji and Si\Ji. (The area Qj is a subset of this one, as known from Lemma^) 
The overcompressive SDSLq, C > 0, is a part of the SDSL bounded by the curves D 
and E such that ui < uq — 2. 

4.1. Interaction of a singular shock and an admissible shock wave. Suppose 
that a singular shock wave with a speed ci and a left- and right-hand values Uo = 
(uo,wo) and C/i = respectively, interact with an admissible shock wave 

with a speed C2 < ci having left-hand and right-hand values Ui — (mi,wi) and 
U2 — (1*2,^2), respectively, at a point (A, T). 

Lemma 2. If the above singular shock and shock wave are admissible, (u2,V2) lies 
between the lines D and E. Thus, the solution after the interaction is a single 
overcompressive singular shock wave. 
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Proof. Since uq > mi + 3 and ui > U2 (because of the admissibility conditions for 
singular and shock wave), we have uq > U2 + 3. The point (u2, V2) lies on the curve 
or S2 with the origin at the point {ui,Vi). Thus 



The point {ui,vi) lies in the area denoted by Qt, thus bellow or at the curve D 
with the origin at {uo,vq). Therefore 



Let the point {uo,vo) be the origin. The point {u2,V2) will be bellow the curve D 



The left-hand side of the above inequality is less than 2, while the right-hand side 

is greater that 2. Thus, the point (u2, ^2) really lies bellow the curve D. 

In the same way one can prove that the point (m2. ^2) lies above the curve E. □ 

Remark 6. In the same manner as above, one can prove that the situation is the 
same when singular shock and shock wave change sides. That is, when an admissible 

singular shock wave interacts with an admissible shock wave from the right-hand 
side, then the solution is again a single admissible singular shock wave. 

4.2. Double singular shock wave interaction. Suppose that an admissible sin- 
gular shock wave with a speed ci and left- and right-hand side values Uq = (uq, vq) 
and Ui = (ui,fi), respectively, interacts with an another singular shock wave with 
a speed C2 < ci and left-hand (right-hand) side values Ui = {ui,vi) {U2 = {u2, V2)) 
at the point {X, T). Since the conditions for the existence of singular shock waves 
include uq — u\> 2> and Ux — U2> 3, then uq — U2 > 6, i.e. the point (u2, V2) is 
on the left-hand side of the line u = uq — \/V2. Concerning the position of the 
point (u2, ^2) in the plane of wave regions with the origin at {uq, vq) there are three 
possibilities: 

(i) The point {u2,V2) is between or at the curves D and E. The result of the 
interaction is a single singular shock wave (with increasing strength). 

(ii) The point (u2,f2) is above the curve D. The result of the interaction is an 
1-rarefaction wave followed with a singular shock wave. 

(iii) The point (u2,f2) is bellow the curve E. The result of the interaction is a 
singular shock wave followed by a 2-rarcfaction wave. 

SDSL's always have increasing strength in these three cases. 

5. Intersection of a singular shock wave and a rarefaction wave 

The last possibility of singular shock wave interactions is between a singular 
shock wave and a rarefaction wave. That possibility is omitted from a considerations 
of the general case due to a richness of possible behaviors. Nevertheless, the most 




Vi <Vo + u{ + Ui- UqUi - Uq. 



if 




<Vo -\-u\+U2 — U0U2 — Uq. 

Non-positivity of Ui — U2 gives 
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of specific Riemann problems can be treated similarly as system ^ was here, at 
least up to some point. 

For a given point (uq, vo)i the rarefaction curves are given by (see (SI) 

Ri = {(u, v) : V = vq- + + u - uq}. 

R2 = {(u, v) : V ^vo- ^Mq + - u + uo} 

Suppose that a singular shock wave with left- and right-hand side values Uq = 
(uo, vq) and Ui = {ui, vi), from the left-hand side interacts with a rarefaction wave 
at some point (X, T) . If the rarefaction wave is approximated with a number of 
small amplitude (non-admissible) shock waves like in wave fronth tracking algorithm 
(see ^ for example), intuition given in Theorem ^ such that the first task should 
be to look at the singular shock wave and the interaction of singular shock and 
non- admissible shock wave. It is possible to extend Theorem ^ for such a case, 
providing that a non-admissible shock wave has amplitude small enough (of the 
rate e^, say). Denote by {ur,Vr) the end-point in a rarefaction curve. Let us note 
that the starting point of the curve {ui,Vi) is in Q7. 

In what follows, we shall abuse the notation and denote by (mi, Vi) the left-hand 
side of an approximated non-admissible shock wave. Denote by {ui,vi) e Q7 the 
left-hand side and by (m2, V2) the right-hand side value of a part from the rarefaction 
curve. If (it2,i'2) G Q7, then the result of the interaction is a single singular shock 
wave, with the left-hand side value equals {uq,vq). The speed depends on initial 
values as in 126|l . So, one can continue the procedure taking approximate points 
from the rarefaction curve as the right-hand values of the non-admissible shock 
wave until it reaches the border of Qy. 

After looking at the above discrete model we are back in a real situation. 

Let us denote by (c(t),t), t belonging to some interval, a path of the resulting 
singular shock wave trough Qy. It is possible to explicitly calculate the above path. 
For example if a singular shock wave interacts with a centered 1-rarefaction waves, 
substituting 

f Ui,a;<c(<)1 
(piix/t), X > c(<) J 

{Vl, X < c(t)l 
^ \ + mDeix - Cit)) 

(p2[x/t), X > C(t) J 

in system one obtains 

c{t) = {t{l - 2{ui -VQ + vi+ul- uf)) 
+ T{1 - 2{uo -Vl- uoui +ul- u?))) /{2{uo - 1)) 
/3'(<) = c\t){l{^ + 1) + (f^ + 1) + _ 



(Kf-r-(f-)-ro-o), 
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where the initial data for j3 at the point t — T \s the initial strength of the singular 
shock wave P{T). The above calculations means that a form of the resulting sin- 
gular shock curve and its strength are uniquely determined trough the area Qr. If 
{ur, Vr) G Qr, then the analysis is finished. Suppose that this is not true. The main 
problem is to analyse situation when rarefaction curve intersects the boundary of 
Qr- Let us try to find out what is happening by using a discrete model. 

Thus, the first real problem is to find a form of solution when the points from 
the rarefaction curve satisfy: e Qr and {u2,V2) ^ Qr- 

Denote by D and G the intersection points of the curve Ji (or the line u = uq — 3) 
with the curves E and D, respectively (see Figure 2). 

5.1. The first critical case. Denote by J the 1-rarefaction curve starting from 
the point G and by J2 the 2-rarefaction curve starting from the point D The region 

where {u2,V2) can lie consist of five subregions: 

(i) The rarefaction curve which starts at {ui,Vi) intersects the curve D out of point 
G. The point (^2,^2) lies in the region above the curve D and left of the line 
w = uo — 3. The final result of the interaction is a 1-rarefaction wave (Ri) followed 
by a singular shock wave with increasing strength. 

(ii) The rarefaction curve which starts at {ui,vi) intersects the curve E out of point 
D. The point {u2,V2) lies in the region below the curve E and on the left-hand side 
of the line u = uq — 3. The result of the interaction is a singular shock wave with 
increasing strength followed by a 2-rarefaction wave (i?2)- 

(iii) The rarefaction curve which starts at (ui,ui) intersects the curve Ji out of 
points D and G. Since an amplitude of a non-admissible shock wave can be as 
small as necessary, one can assume that the point {u2, V2) lies in the second delta 
singular locus and the resulting singular shock wave has a negative strength. The 
strength-function /?(t) = + k{t — To) of the resulting singular shock is decreasing, 
so, there could exists a point Ti — T — C^/k such that 0{Ti) = ao = ai = 0. Let 
Xi = cTi -|- {X — T), where c is the speed of the resulting singular shock wave 
(space coordinate of the point where strength reaches zero). Therefore, in the time 
t — Ti, we have to solve new Riemann problem 



This problem has a unique entropy solution consists from two shock waves, since 
the point (1*2,^2) is between the curves 5*1 and 5*2, with respect to the origin at 
the point {uq,vq). This means that the singular shock wave decouples into a pair 
of admissible shock waves. If <uq—2, this pair of the shock waves are the final 
solution. The case when Ur < uq — 2 belongs to the following subsection, i.e. the 
second critical case. 

(iv) The rarefaction curve Rj, j = 1 or 2, which starts at {ui,vi) intersects the 
curve Ji in the point G- We can take G = {U21V2) for convenience. The set of such 
points (ui, Ul) lies on the inverse rarefaction curve, which starts from the right-hand 
side values, i.e. 




X < Xi 
x> Xi 



Ri = {{u, v) V = vo + (ul — M^)/2 + uo — u} 
and 



R2 = {{u,v) : V = Vo + {ul - u^)/2 -uo + u} 
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(the same explanation will be used in Remark|Z|bcllow) . Straightforward calculation 
shows that this curve lies in the region Qj, thus this situation is possible, as one 
can see using the inverse rarefaction curves Ri and R2 given above. 

If J — 1, then the point [u2, V2) belongs to J and the solution after the interaction 
is an i?i-wave followed by a singular shock wave with a constant strength. 

If j = 2, then the point (^2,^2) lies in the area bellow the curve J. This can 
be verified by direct calculation, taking into account that the amplitude of a non- 
admissible shock is small enough, U2 < uo — 2. The solution after the interaction is 
an admissible singular shock wave with a decreasing strength. Further explanations 
of a such singular shock wave is given in the following subsection, 
(v) The rarefaction curve Rj which starts at (ui,vi) intersects the curve Ji in the 
point D. Again, let D = {u2,V2)- Simple calculation, as in the case (iv), shows 
that this situation is also possible since the inverse rarefaction curves starting from 
D stay in Q-;. If j — 2, the point (u2,i'2) belongs to J2, and then the solution 
after the interaction is singular shock wave with a constant strength followed by an 
i?2-wave. 

If j = 1, then use of the same arguments as above gives that the point (u2, "^2) 
lies in the area above the curve J2 and the result of the interaction is an admissible 
singular shock wave with a decreasing strength. Again, one can see the following 
subsection for the further analysis. 

5.2. The second critical case. Now we are dealing with the problem when the 
rarefaction wave after passing through Ji, after passes trough the curves D or E. 
One can see that this is the continuation of the cases (iii)-(v) from the previous 
part. 

(a) Denote by D the area above the curve Z?, bellow 5*1 and on the left-hand side of 

the line u = uq — 2. Also denote by D the area above the curve E, bellow Si and on 
the right-hand side of the line u = uq — 2. If (u2, ^2) lies in one of these regions, the 
solution is combination of a rarefaction wave i?i and an overcompressive singular 
shock wave with a decreasing or constant strength. 

(b) Denote by E the area bellow the curve E^ above S2 and on the left-hand side 

of the line u ~ — Also denote by E the area bellow the curve Z), above S'2 and 
on the right-hand side of the line u = — 2. If (7/2, ^2) lies in one of these regions, 
the solution is then a combination of an overcompressive singular shock wave with 
a decreasing or constant strength and a rarefaction wave i?2. Denote by Dq the 
area bounded by the curves D, E, Si and S'2 such that u < wq — 2 in Dq- 

One can see that a rarefaction curve cannot enter into Dq since it has to pass 
trough the intersection point {uq — 2, —2uq + vq + 2) of I? and E, but 

2u — 2uo — uuq + u'^/2 + Uq/2 + 2 > (ie. Ri is above the curve E) 
2uo — 2u + uuo — V? 12 — Uq/2 — 2 < (ie. R2 is bellow the curve D). 

Therefore, a rarefaction curve which passes trough the point DClE goes either into 

D or E, and these cases are analysed above. 

Thus, we have described all important points of the interactions between singular 
shock and rarefaction waves. When a result of a single interaction is known, the 
question about further singular shock path could be answered by a successive use 
of the above procedures. 
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Figure 2. 



Remark 7. One can use the similar analysis of all possible cases when a rarefaction 
wave which interacts with a singular shock wave is on the left-hand side of it. 
Instead of direct rarefaction and singular shock curves, the inverse ones should be 
used, i.e. (it2,i'2) is a starting point and one is able to calculate vq from formulas 
ofE ,D, Si, 52, Ri and i?2- 

Remark 8. In the contrast with the case in 0, where interaction can generate some 
"strange" solution containing unbounded L}^^ function, in the presented system 
one can find only bounded functions and singular shock waves as a result on an 
interaction. 

For a system ^ with gi ^ const, or 52 ^ 0; interaction of singular shock and 
rarefaction waves cannot be treated as easy as here. 

Thus, we have proved the following assertion for the interaction in the case of 
system 

Theorem 3. Suppose that a singular shock wave interacts with a rarefaction wave 
at the time T . For some time period T < t < Ti the solution is represented by a 
singular shock wave supported by a uniquely defined curve (not a line) followed by a 
new rarefaction wave. Depending on the right-hand value of the primary rarefaction 
wave, one has the following possible cases for a solution after t > Ti. 

(a) Single singular shock wave ( supported by a line ) with an increasing strength. 

(b) 1-rarefaction wave followed by singular shock wave with an increasing strength. 

(c) Singular shock wave with an increasing strength followed by 2-rarefaction 
wave. 
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(d) Singular shock wave with a decreasing strength prolonged by either a single 
singular shock wave with an increasing strength, or a pair of admissible 
shock waves. 

(e) 1-rarefaction wave followed by singular shock wave with a constant strength. 

(f) Singular shock wave with a constant strength followed by 2-rarefaction wave. 

(g) Singular shock wave with a decreasing strength prolonged by either 1-rarefaction 
wave followed by singular shock wave with decreasing or constant strength, 

or singular shock wave with decreasing or constant strength followed by 2- 
rarefaction wave. 

"Prolonged" is the state after strength of singular shock wave becomes zero. Such 
wave can also stop with non-zero strength, and then there is obviously no prolonga- 
tion described above. 
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